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ABSTRACT 

We study the density field around Zem > 4 quasars using high quality medium spectral reso- 
lution ESI-Keck spectra {R ^ 4300, SNR > 25) of 45 high-redshift quasars selected from a 
total of 95 spectra. This large sample considerably increases the statistics compared to previ- 
ous studies. The redshift evolution of the mean photo-ionization rate and the median optical 
depth of the intergalactic medium (IGM) are derived statistically from the observed transmit- 
ted flux and the pixel optical depth probability distribution function respectively. This is used 
to study the so-called proximity effect, that is, the observed decrease of the median optical 
depth of the IGM in the vicinity of the quasar caused by enhanced photo-ionization rate due 
to photons emitted by the quasar We show that the proximity effect is correlated with the 
luminosity of the quasars, as expected. By comparing the observed decrease of the median 
optical depth with the theoretical expectation we find that the optical depth does not decrease 
as rapidly as expected when approaching the quasar if the gas in its vicinity is part of the 
standard IGM. We interpret this effect as revealing gaseous overdensities on scales as large 
as ~I5/i~^ Mpc. The mean overdensity is of the order of two and five within, respectively, 
10 and 3/i~^ Mpc. If true, this would indicate that high redshift quasars are located in the 
center of overdense regions that could evolve with time into massive clusters of galaxies. The 
overdensity is correlated with luminosity: brighter quasars show higher overdensities. 

Key words: Methods: data analysis - statistical - Galaxies: clustering - intergalactic medium 
- quasai^s: absorption lines - Cosmology: dai^k matter 



1 INTRODUCTION 

The Inter-Galactic Medium (IGM) has been intensively studied 
using the absorption seen in the spectra of quasi-stellar objects 
(QSOs) over a large redshift range (0.16 < Zcm < 6.3). This ab- 
sorption, first identified by Lynds (1971), breaks up at high spectral 
resolution in hundreds of discrete absorption lines from, predom- 
inantly, HI Lyman UV resonance lines redshifted in a expanding 
universe (the so-called Lya forest, see Ranch 1998 for a review). 

The Lya forest was interpreted by Sargent et al. (1980) as the 
signature of intervening HI clouds of cosmological nature embed- 
ded in a diffuse hot medium. The clouds were further described 
by Rees (1986) as gravitationally confined by dark-matter mini- 

* Based on observations carried out at the Keck Telescope 



halos. The advent of numerical simulations has introduced a new 
and more general scheme in which the IGM is a crucial element of 
large scale structures and galaxy formation. It is now believed that 
the space distribution of the gas traces the potential wells of the 
dark matter. In addition, most of the baryons are in the IGM at high 
redshift, making the IGM the reservoir of gas for galaxy formation. 
The numerical A'^-body simulations have been successful at repro- 
ducing the observed characteristics of the Lya forest (e.g., Cen et 
al. 1994; Petitjean et al. 1995; Hemquist et al. 1996; Theuns et al. 
1998). The IGM is therefore seen as a smooth pervasive medium 
which can be used to study the spatial distribution of the mass on 
scales larger than the leans' length. This idea is reinforced by ob- 
servations of multiple lines of sight (e.g., Coppolani et al. 2006). 

It is well known that the characteristics of the Lya forest 
change in the vicinity of the quasar due to the additional ionizing 
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flux produced by the quasar. The mean neutral hydrogen fraction 
decreases when approaching the quasar. Because the amount of ab- 
sorption in the IGM is, in general, increasing with redshift, this re- 
versal in the cosmological trend for redshifts close to the emission 
redshift of the quasar is called the 'inverse' or 'proximity' effect 
(Carswell et al. 1982; Murdoch et al. 1986). It is possible to use 
this effect, together with a knowledge of the quasar luminosity and 
its position, to derive the mean flux of the UV background if one 
assumes that the redshift evolution of the density field can be ex- 
trapolated from far away to close to the quasar. Indeed, the strength 
of the effect depends on the ratio of the ionization rates from the 
quasar and the UV background, and because the quasar's ionization 
rate can be determined directly through the knowledge of its lumi- 
nosity and distance, the ionization rate in the IGM can be inferred. 
This method was pioneered by Bajtlik, Duncan & Ostriker (1988) 
but more recent data have yielded a wide variety of estimates (Lu, 
Wolfe & Tumshek 1991; Kulkami & Fall 1993; Bechtold 1994; 
Cristiani et al. 1995; Femandez-Soto et al. 1995; Giallongo et al. 
1996; Srianand & Khare 1996; Cooke, Espey & Carswell 1997; 
Scott et al. 2000, 2002; Liske & Williger 2001). Scott et al. (2000) 
collected estimates from the literature which vary over almost an 
order of magnitude at z = 3. The large scatter in the results can be 
explained by errors in the continuum placement, cosmic variance, 
redshift determination, etc. 

In the standard analysis of the proximity effect it is assumed 
that the underlying matter distribution is not altered by the pres- 
ence of the quasar. The only difference between the gas close to the 
quasar or far away from it is the increased photoionization rate in 
the vicinity of the QSO. If true, the strength of the proximity effect 
should correlate with the quasar luminosity but such a correlation 
has not been convincingly established (see Lu et al. 1991; Bech- 
told 1994; Srianand & Khare 1996; see however Liske & Williger 
2001). It is in fact likely that the quasars are located inside over- 
dense regions. Indeed, the presence of Lya absorption lines with 
2abs > «em suggests an excess of material around QSOs (Loeb & 
Eisenstein 1995; Srianand & Khare 1996). Furthermore, in hierar- 
chical models of galaxy formation, the supermassive black holes 
that are thought to power quasars are located in massive haloes 
(Magorrian et al. 1998; Ferrarese 2002), that are strongly biased 
to high-density regions. Possible evidence for overdensities around 
quasars come also from studies of the transverse proximity effect 
by Croft (2004), Schirber, Miralda-Escude & McDonald (2004) 
and Worseck & Wisotzki (2006) who suggest that the observed 
absorption is larger than that predicted by models asstuning stan- 
dard proximity effect and isotropic quasar emission. However, in 
the case of transverse observations, it could be that the quasar light 
is strongly beamed in our direction or, alternatively, that the quasar 
is highly variable. Interestingly, neither of these affects the longitu- 
dinal proximity effect discussed in the present paper. 

Observations of the IGM transmission close to Lyman break 
galaxies (LBGs) seem to show that, close to the galaxy, the IGM 
contains more neutral hydrogen than on average (Adelberger et al. 
2003). As the UV photons from the LBGs cannot alter the ion- 
ization state of the gas at large distances, it is most likely that the 
excess absorption is caused by the enhanced IGM density around 
LBGs. It is worth noting however that various hydrodynamical sim- 
ulations have trouble reproducing this so-called galaxy proximity 
effect (e.g., Bruscoh et al. 2003; KoUmeier et al. 2003; MaseUi et 
al. 2004; Desjacques et al. 2004). 

In a recent paper, RoUinde et al. (2005) presented a new anal- 
ysis to infer the density structure around quasars. The method is 
based on the determination of the cumulative probability distribu- 



tion function (CPDF) of pixel optical depth, and so avoids the Voigt 
profile fitting and line counting which is traditionally used (e.g., 
Cowie & Songaila 1998; ElUson et al. 2000; Aguirre et al. 2002; 
Schaye et al 2003; Aracil et al. 2004; Fieri et al. 2006). The evolu- 
tion in redshift of the optical depth CPDF far away from the quasar 
is directly derived from the data. This redshift dependent CPDF is 
then compared to the CPDF observed close to the quasar to derive 
the mean density profile around quasars. The method was applied 
to twenty lines of sight toward quasars at Zam ~ 2 observed with 
UVES/VLT and it was found that overdensities of the order of a 
few are needed for the observations to be consistent with the value 
of the UV background flux derived from the mean Lya opacity. In 
the present paper we apply the same method to a large sample of 
95 quasars at Zom > 4 observed with the Echelle Spectrograph and 
Imager (ESI) mounted on the Keck II telescope. 

In Section 2 we describe the data and the selection of the sam- 
ple used in the present work. We derive the redshift evolution of 
the ionizing UV background and of the median IGM optical depth 
in Sections 3 and 4 respectively. We discuss the proximity effect 
in Section 5 and conclude in Section 6. We assume throughout this 
paper a flat Universe with fim = 0.3, Ha = 0.7, Ob = 0.04 and 
H<, = 70kms-^ 



2 DATA AND SAMPLE SELECTION 

Medium resolution (R ~ 4300) spectra of all z > 3 quasars 
discovered in the course of the DPOSS survey (Digital Palomar 
Observatory Sky Survey; see, e.g., Kennefick et al. 1995, Djor- 
govski et al. 1999 and the complete listing of QSOs available at 
http://www.astro.caltech.edu/'N^george/z4.qsos) have been obtained 
with the Echellette Spectrograph and Imager (ESI, Sheinis et al 
2002) mounted on the KECK II 10 m telescope. Signal-to-noise ra- 
tio is usually larger than 15 per 10 km s"'^ pixel . These data have 
already been used to construct a sample of Damped Lya systems 
at high redshift (Prochaska et al. 2003a,b). In total, 95 quasars have 
been observed. 

The KECK/ESI spectra were reduced (bias subtraction, flat- 
fielding, spectrum extraction) using standard procedures of the 
IRAF package. The different orders of the spectra were combined 
using the scombine task. In the Echellette mode spectra are di- 
vided in ten orders covering the wavelength range: 4000 A < 
Aobs < 10000 A. In this work we use only orders 3 (center 4650 A) 
to 7 (center 6750 A). The spectral resolution is i? ^ 4300 or 
~ 70 km s^^. During the process of combining the orders, we 
controlled carefully the signal-to-noise ratio obtained in each or- 
der. Wavelengths and redshifts were computed in the heliocentric 
restframe and the spectra were flux calibrated and then normalized. 

The signal-to-noise ratio per pixel was obtained in the regions 
of the Ly-a forest that are free of absorption and the mean SNR 
value, averaged between the Ly-a and Ly-/3 emission lines, was 
computed. We used only spectra with mean SNR > 25. We rejected 
the broad absorption line QSOs (BAL) and QSOs with more than 
one damped Lya system (DLA) redshifted between the Ly-a and 
the Ly-/3 emission lines because their presence may over-pollute 
the Lya forest. Metal absorptions are not subtracted from the spec- 
tra. We can estimate that the number of intervening CIV and Mgll 
systems with Wobs > 0.25 A is of the order of five along each 
of the lines of sight (see e.g., Boksenberg et al. 2003; Aracil et al. 
2004; Tytler et al. 2004; Scannapieco et al. 2006). This means we 
expect an error on the determination of the mean absorption of the 
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Figure 1. Number of spectra in our sample contributing to the study of the 
Ly-a forest at a given redshift as a function of redsliift. 



order of 1 %. This is at least five times smaller than the error ex- 
pected from the placement of the continuum. 

In Table 1 we give the characteristics of the forty-five QSO 
spectra satisfying the above criteria that are used in the present 
work. Column (1) gives the QSO's name; column (2) the emis- 
sion redshift estimated as the average of the determinations of the 
peak of the Lya emission and the peak of a Gaussian fitted to the 
CIV emission line; column (3) the emission redshift obtained using 
the IRAF task rvidlines (see Section 5.1); column (4) the appar- 
ent V magnitude; column (5) the mean signal-to-noise ratio in the 
Lya forest and column (6) the intrinsic luminosity at the Lyman 
limit estimated from the V-magnitude assuming that the QSO con- 
tinuum spectrum is a power-law of index —0.6 (see Section 5). In 
Fig. 1 we show the histogram of the number of spectra in our sam- 
ple contributing to the study of the Lya forest at a given redshift as 
a function of redshift. Normalisation of the spectra is known to be a 
crucial step in these studies. An automatic procedure (Aracil et al. 
2004) estimates itcrativcly the continuum by minimising the sum 
of a regularisation term (the effect of which is to smooth the con- 
tinuvim) and a term, which is computed from the difference be- 
tween the quasar spectrum and the continuum estimated during the 
previous iteration. Absorption lines are avoided when computing 
the continuum. A few obvious defects are then corrected by hand 
adjusting the reference points of the fit. This happens to be impor- 
tant near the peak of strong emission lines and over damped ab- 
sorption lines. The automatic method works very well and a mini- 
mal manual intervention is necessary. The procedure was calibrated 
by Aracil et al. (2004) using simulated quasar spectra (with emis- 
sion and absorption lines) adding continuum modulations to miiruc 
an imperfect correction of the blaze along the orders and noise to 
obtain a S/N ratio similar to that in the data. We noted that the pro- 
cedure underestimates the true continuum in the Lya forest by a 
quantity depending smoothly on the wavelength and the emission 
redshift by an amount of about 3-5% at 2: ~ 3.5-4. (see Fig. 1 of 
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Aracil et al. 2004). This is less than our typical errors and therefore 
we did not correct the normalized spectra for this. Note that, due 
to strong blending, errors can be as high as 10% in places (see also 
Croft et al 2002; Becker et al. 2006; McDonald et al. 2006; Des- 
jacques, Nusser & Shcth 2007). In that case however, the optical 
depth will be usually larger than the limit we will used to reject 
pixels badly affected by saturation. It must be noted also that, be- 
cause of the reduced absorption in the vicinity of the quasar, the 
continuum determination is more reliable in the wing of the Lya 
emission line. The normalization procedure we have used is pre- 
ferred to other methods that are more arbitrary. It is completely 
reproducible and introduces only a small systematic error that we 
can control. A final sample reduced spectrum with the continuum 
fitting is shown in Fig. 2. 



3 THE IONIZING BACKGROUND RADIATION FIELD 

The standard treatment of the proximity effect consists in studying 
the evolution of the mean absorption in the Lya forest when ap- 
proaching the quasar. Far away from the quasar the only source of 
ionizing photons is the UV background whereas in the vicinity of 
the quasar, the gas is ionized by both the UV background and the 
QSO. Assuming (i) that the luminosity of the quasar is known, (ii) 
that the distance from the gas to the quasar is cosmological (and 
therefore given by the difference in redshift) and (iii) that the den- 
sity field in the IGM is not modified by the presence of the quasar, 
it is possible to derive observationally the distance to the quasar 
where the ionizing flux from the quasar equals the flux from the UV 
background. This, in turn, gives an estimate of the UV backgroimd 
flux. The last assumption neglects the fact that quasars can be sur- 
rounded by significant overdensities (e.g., PascareUe et al. 2001; 
Adelberger et al. 2003; RoUinde et al. 2005; Faucher-Giguere et al. 
2006; Kim & Croft 2006). 

However the above approach can be reversed to derive the den- 
sity distribution around the quasar if the UV background can be 
estimated from elsewhere. Actually, it is possible to estimate this 
flux by modelling the redshift evolution of the mean absorption in 
the IGM. For « > 4 we can follow Songaila & Cowie (2002) and 
use the mean normalized transmitted flux, F{z), to derive the nor- 
malized ionization rate, g, in units of I0~^^ s~^, by inverting the 
following equation (see also McDonald & Miralda-Escud6 2001; 
Cen & McDonald 2002): 

F{z) = 4.5 X X [^^r X el--^^-""!^!^! (1) 

with 

g = r_i2 X t!-'' xhx X (2) 

^ * ^0.35^ ^0.0325^ 

Applying the above equations, we can derive the redshift de- 
pendence of r_i2 from the observed transmitted flux (for z > 4). 
For this, we have used the Lya forest over the rest-wavelength 
range 1070-1170 A to avoid contamination by the proximity ef- 
fect close to the Lya emission line and by possible OVI associated 
absorbers. We also carefully avoided regions flagged because of 
data reduction problems or damped Lya absorption lines. We di- 
vided each spectrum in bins of length 50 A in the observed frame, 
corresponding to about Az = 0.04. The mean transmitted flux 
was calculated as the mean flux over all pixels in a bin. At each 
redshift, we then averaged the transmitted fluxes over all spectra 
covering this redshift. Errors were estimated as the standard devi- 
ation of the mean values divided by the square root of the number 
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Figure 2. Flux calibrated spectrum of PSS 1058+1245. The continuum derived automatically is overplotted. 
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Figure 3. The observed normalized transmitted flux < F > as a function 
of redshift. 



of spectra. The corresponding scatter in the transmitted flux cannot 
be explained by photon noise or by uncertainties in the continuum, 
suggesting that errors are dominated by cosmic variance. The mean 
transmitted flux that we have obtained is plotted in Fig 3 together 
with results by Songaila (2004). 

The normalized ionization rate defined in Eq. 2 is derived from 
the mean normalized transmitted flux by inverting Eq. 1. Results 
are plotted in Fig. 4. To decrease the errors, we have used only three 
bins at 2 > 4. It can be seen that our measurements are consistent 
with those by McDonald et al. (2000) at lower redshift. Fitting both 
McDonald's and our results together, we find that the redshift evo- 
lution of g is described as 3 = l.llx[(l+z)/6]~^'^^. Note that our 
fit is not consistent with the results by Songaila (2004) who find 
g = 0.74x [(l+z)/6]^^'^ for z > A. The corresponding photoion- 
ization rate, r_i2, in units of 10^^'^ s^^, is given in Fig. 5, assum- 
ing Q,^h^ = 0.019, fim = 0.3 and three different values of gas 
temperature r4 = 1, 1.5 and 2, in units of 10* K. These results are 
consistent with the measurements by McDonald & Mirada-Escude 
(2001) for a mean temperature of r4 1.5 to 2 which is the mean 
temperature expected in the IGM at these redshifts. However, re- 
cent determination of this quantity (Becker et al. 2006) using a log- 



Figure 4. The normalized ionization rate obtained from the transmitted flux 
observed in our data by inverting Eq. 1 in three redshift bins (filled squares) 
at > 4. Filled circles are the measiu'ements of McDonald et al. (2000) 
for lower redshift. The dotted line is a power law of the form 1.11 x [( 1 + 
z) /6] -1-83 flttgd to flie data of McDonald et al. (2000) and fliis paper. The 
dashed line is a power law of the form 0.74 x [(1 -|- z)/&\~'^'^ obtained, 
for z > 4, by Songaila (2004). 



normal distribution for the optical depth distribution indicate that 

r-12 could be higher by about a factor of two to three at these red- 
shifts. We therefore have used in the following a mean temperature 
0fT4 = 1. 



4 THE Lya OPTICAL DEPTH STATISTICS IN THE IGM 

To study the influence of the additional ionizing flux from the 
quasar on the Lya optical depth, we have first to derive the evolu- 
tion of the optical depth with redshift in the IGM at large distances 
from the quasar. The observed HI opacity is simply: 



rHi(A) = -ln( 



Fx.. 



(3) 



is the flux in the con- 
tinuum. Its evolution with redshift is described as 



where ^x^bs observed flux and J^Aeont 
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Table 1. Quasar sample 



QSO 




2^em 


Vmag 


SNR 


^912 

^, — 2 

n ergs s 




PSSOl 17+1552 


4.241 


4.244 


18.6 


45 


4.53 X 


1031 


PSS0118+0320 


4.235 


4.232 


18.50 


30 


4.38X 


1031 


PSS0121+0347 


4.130 


4.127 


17.86 


72 


6.76X 


1031 


SDSSO 127-0045 


4.067 


4.084 


18.37 


25 


4.04 X 


1031 


PSSOl 3 1+0633 


4.432 


4.430 


18.24 


25 


6.88X 


1031 


PSSOl 34+3307 


4.534 


4.536 


18.82 


30 


4.73 X 


1031 


PSS0209+0517 


4.206 


4.194 


17.36 


40 


8.16X 


1031 


PSS0211+1107 


3.975 


3.975 


18.12 


66 


6.14X 


1031 


PSS0248+1802 


4.427 


4.430 


18.4 


60 


8.32X 


1031 


PSS0452+0355 


4.397 


4.395 


18.80 


30 


4.77 X 


1031 


PSS0747+4434 


4.434 


4.435 


18.06 


54 


7.98X 


1031 


PSS0808+5215 


4.476 


4.510 


18.82 


37 


4.29 X 


1031 


SDSS0810+4603 


4.078 


4.074 


18.67 


43 


3.22X 


1031 


PSS0852+5045 


4.213 


4.216 


19 


54 


2.60 X 


1031 


PSS0926+3055 


4.188 


4.198 


17.31 


77 


1.22X 


1032 


PSS0950+5801 


3.969 


3.973 


17.38 


85 


9.31 X 


1031 


PSS0957+3308 


4.283 


4.274 


17.59 


40 


9.80 X 


1031 


PSS 1057+4555 


4.127 


4.126 


17.7 


76 


7.67 X 


1031 


PSS1058+1245 


4.332 


4.330 


18 


46 


7.45 X 


1031 


PSS 1140+6205 


4.507 


4.509 


18.73 


61 


4.77 X 


1031 


PSS1159+1337 


4.089 


4.081 


18.50 


45 


3.65X 


1031 


PSS 1248+3 110 


4.358 


4.346 


18.9 


25 


3.19X 


1031 


SDSS1310-0055 


4.152 


4.151 


18.85 


38 


2.82X 


1031 


PSS1317+3531 


4.370 


4.369 


19.10 


28 


2.75X 


1031 


J1325+1123 


4.408 


4.400 


18.77 


25 


4.20 X 


1031 


PSS1326+0743 


4.121 


4.123 


17.3 


66 


1.01 X 


1032 


PSS1347+4956 


4.597 


4.560 


17.9 


40 


1.02 X 


1032 


PSS1401+4111 


4.008 


4.026 


18.62 


30 


2.90 X 


1031 


PSS1403+4126 


3.862 


3.862 


18.92 


25 


1.85 X 


1031 


PSS 1430+2828 


4.309 


4.306 


19.30 


40 


2.18X 


1031 


PSS 1432+3940 


4.291 


4.292 


18.6 


36 


4.04 X 


1031 


PSS 1443+2724 


4.419 


4.406 


19.3 


30 


2.57 X 


1031 


PSS 1458+68 13 


4.295 


4.291 


18.67 


60 


5.04 X 


1031 


PSS 1500+5829 


4.229 


4.224 


18.6 


40 


3.74 X 


1031 


GB1508+5714 


4.306 


4.304 


18.9 


32 


3.12X 


1031 


PSS1535+2943 


3.979 


3.972 


18.9 


29 


2.26 X 


1031 


PSS1555+2003 


4.226 


4.228 


18.9 


31 


3.15X 


1031 


PSS1633+1411 


4.360 


4.349 


19.0 


43 


3.36 X 


1031 


PSS1646+5514 


4.110 


4.084 


18.11 


37 


4.81 X 


1031 


PSS1721+3256 


4.031 


4.040 


19.23 


43 


1.85X 


1031 


PSS 1723+2243 


4.515 


4.514 


18.17 


42 


8.77 X 


1031 


PSS2154+0335 


4.349 


4.359 


18.41 


26 


6.20 X 


1031 


PSS2203+1824 


4.372 


4.375 


18.74 


34 


4.50X 


1031 


PSS2238+2603 


4.023 


4.031 


18.85 


24 


2.74X 


1031 


PSS2344+0342 


4.341 


4.340 


17.87 


30 


4.38 X 


1031 



" Mean of estimates from a Gaussian fit to the CIV emission line and from 
the peak of the Lya emission 
^ Estimate using the IRAF RVIDLINES task 
continuimi luminosity at 912 A 



Toc(l + z)°. (4) 

For each spectrum of the sample we estimate the optical depth, 
T, in each pixel between Tmin = — ln(l — 3a) and Tmax = 
— ln(3(T), where tj(A) is the rms of the noise measured in the spec- 
trum (see column 5 of Table 1). We then construct the cumulative 
probability distribution of r (CPDF) in redshift windows of approx- 
imately Az = 0.04 corresponding to 50 A in the observed frame 



and estimate the associated percentiles. Because low and high val- 
ues of the optical depth are lost either in the noise or because of sat- 
uration, we can use only the intermediate values of the percentiles. 
The evolutions with redshift of the 40, 50, 60 and 70 % percentiles 
are given in Fig. 6. 

The redshift evolution index in Eq. 4 can be derived from the 
redshift evolution of the percentiles of the pixel optical depth CPDF 
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Lognormal model 
McDonald & Miralda-Escude (2001) 
This work T=1 OOOOK 
This work T=1 5000K 
This work T=20000K i — 



40th percentile I \ — 

50th percentile ! — x— 

60th percentile 

70th percentile r B -- 



Figure 5. Redshift evolution of the photoinization rate, r_ 12, derived from 
our data and using three different temperatures for the IGM, T4 = 1.0, 
1.5 and 2, in units of 10^ K (filled diamond, open diamond and crosses, 
respectively). Filled squares show the results by McDonald & Miralda- 
Escude (2001) and filled circles the results from Becker et al. (2006). 
Q,\,H^ = 0.019, Qm = 0.3 are assumed here. 



(see RoUinde et al. 2005). The values obtained for a are (5.8 ± 0.5) 
for the 40 % percentile, (5.1 ± 0.2) for the 50 % percentile,(4.3 ± 
0.2) for the 60 % percentile and (3.8 ± 0.4) for the 70 % percentile. 
We use throughout the rest of the paper a=4.8. 



5 THE PROXIMITY EFFECT FROM OPTICAL DEPTH 
STATISTICS 

The evolution with redshift of the optical depth in the IGM is in- 
verted in the vicinity of the quasar due to the additional ioirizing 
flux from the quasar. To disentangle the two effects we can correct 
the observed optical depth for the redshift evolution in the absence 
of the quasar by replacing t\ at redshift z\ by n x [ ] " , where 

Zjci is a fixed redshift taken as reference. We have corrected the 
observed optical depth for each pixel using the factor above with 
Ziei = 4.5 and the value of a obtained from the redshift evolution 
of the IGM CPDF, a = 4.8. The results are shown in Fig 7 (filled 
squares). The observed corrected median optical depth is given ver- 
sus the distance to the quasar computed using the cosmological pa- 
rameters given in Section 2 and assuming emission redshifts from 
column 3 of Table 1. The proximity effect is apparent as a decrease 
of the optical depth within a distance to the quasar smaller than 
~15-20 /i-^Mpc. 



Figure 6. Redshift evolution of different percentiles derived from the CPDF 
of the pixel optical depth. 
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Figure 7. The median optical depth corrected for redshift evolution (t cx 

(1 + 2)" with a = 4.8) is given versus the distance to the quasar, derived 
using QSO emission redshifts from column 3 of Table 1 (filled squares). 
Open triangles correspond to the case where these redshifts are increased 
by an amount taken randomly between and 1500 km s~^. 



5.1 QSO systemic redshifts 

The above effect depends on the accurate determination of the QSO 
systemic redshift. Gaskell (1982) has shown that the redshifts de- 
rived from different quasar emission lines often do not agree with 
each other within typical measurement errors. The high-ionization 
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broad emission lines (HILs; e.g., Lya A1215.67, CIV A1549, CIII 
A1909 and NV A1240) are found systematically blueshifted with 
respect to the low-ionization broad emission lines (LILs; e.g, 01 
A1305, Mgll A2798, and the permitted Balmer series) and princi- 
pally from the forbidden narrow emission lines (e.g., NeV A3426, 
on A3727, Nelll A3870, OHI A5007). The narrow emission lines 
arise from gas located in the galaxy host and their redshift should be 
more representative of the center-of-mass redshift (see e.g., Tytler 
& Fan 1992, Baker et al. 1994). 

Unfortunately, for our sample of quasars with emission red- 
shifts in the range 4 < «em < 4.5, the forbidden narrow emis- 
sion lines cannot be detected from the ground because they are 
redshifted in near-infrared spectral windows that are absorbed by 
the terrestrial atmosphere. We have derived the emission redshifts 
using two approachs. In the first approach, the final redshift (see 
column 2 of Table I) is the mean of two estimates: one obtained 
by fitting a Gaussian to the CIV emission line and the other by 
measuring the peak of the Lya emission line to avoid the various 
absorption features shortwards of the line. In the second approach, 
the emission redshift (see column 3 of Table 1) was obtained using 
the IRAF task rvidlines. Initially we identify a prominent spectral 
feature (usually the Lya emission line) to which a gaussian func- 
tion is fitted. Based on the central wavelength of this line and an 
input list of known spectral features, other features at a consistent 
redshift are identified and fitted. The final emission redshift is a 
weighted average value based on the gaussian fits. The median op- 
tical depth of the IGM versus the distance to the quasar computed 
using emission redshifts obtained by this second approach is shown 
as squares in Fig. 7. 

As we know that HILs are systematically blueshifted, with re- 
spect to LILs and narrow emission lines, from 500 to 1500 km s~^ 
(Tytler & Fan 1992) we have increased the emission redshift (ob- 
tained by fitting the emission lines) by a random amount and cal- 
culated the optical depth versus the distance to the quasar as fol- 
lows. For each realization, we increase the redshift of each of the 45 
quasars by an amount taken randomly between and 1500 km s~^. 
We calculate the distance of each pixel of the lines of sight to the 
corresponding quasar and compute the median optical depth for 
each value of r by averaging over all quasars. We then average 
the optical depths over a hundred realizations. Errors are taken as 
the mean rms obtained over the realizations. Results are shown in 
Fig. 7 (triangles). The proximity effect is more pronounced in that 
case, as expected. This approach will be used in the next Sections 
whenever we will use QSO emission redshifts. 
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Figure 8. The ratio uj = of the ionizing rate from the QSO to 

that from the UV background versus the cosmological distance between 
the quasar and the cloud in the IGM. Tbck is taken from Fig. 5. 



with Lqso the monochromatic limrinosity of the quasar at the Ly- 
man limit. These luminosities are computed extrapolating the flux 
in the continuum at Aobs ~ 6000 A using a power-law of index 
(j) = —0.6 (Francis et al. 1993). We checked that within a reason- 
able range of (j> = -0.5 to -0.7 (e.g., Cristiani & Vio 1990), our 
main result (i.e. the density structure aroimd quasars) is not affected 
by this choice. Therefore, 



Fqso 
thus, 

1^12 
QSO 



12.6 X 10" 

3TT 



-Jqso(i'o) 10^^(s ^) 



12.6 LQso/47r 1^21 (^^-i^^- 



3-1- 



47r 



(7) 



(8) 



r is the luminosity distance from the quasar of emission redshift 
^em to the cloud at redshift z and is calculated using equations 
given by Liske (2003) for a flat cosmological model. 



5.2 The QSO ionization rate 

The strength of the proximity effect depends on the ratio of the ion- 
ization rates from the QSO emission and the UV-background. The 
QSO ionization rate can be determined directly from its luminosity. 
The HI ionization rate due to a source of UV photons is formally 
given by the equation: 



Fqso 



-f 



[47rJQso(i^) 



hv 



(5) 



the frequency of the Lyman limit, CTHi(i- 



cm is the HI photo-ionization cross-section. 



where, vo is 
6.3 X 10-i*(^) 
jQSoiy) = Jq,so{vo) X {-^)~'^, if we assume that the ionizing 
spectrum is a power law of index cj) and where, Jqso {vo) is defined 
as: 

^QSO 



47r jQso(fc'o) 



47^7-2 



(6) 



5.3 Ionizing rate ratio 

We define w as the ratio of the ionizing rate due to the quasar emis- 
sion to that due to the UV background, lj = ""qso ^ j ^ 
has been derived in Section 4 (see Fig. 5). The enhanced ionizing 
flux in the vicinity of the QSO induces a decrease in the optical 
depth of the IGM observed in the absence of the QSO by a factor 
(1 + uj). For each spectrum, we have calculated loqso and r for 
each pixel and then, at a given r, we have derived ui for the whole 
sample as the median of the individual values a;Qso(r) foimd for 
each of the 45 quasars. The result is shown in Fig. 8. 



5.4 Overdensities around the quasar 

In the standard analysis of the proximity effect it is assumed that 
the matter distribution in the IGM is not altered by the presence of 
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density prolil 6 this work 
underestimated redshift and greater ionization rate <, 
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Figure 9. Comparison of median optical depth's evolution observed close 
to the quasars (filled circles) with that of an evolution due to the ionization 
field of the quasar only (Eq. 10; filled squares). The difference between the 
curves can be interpreted as the presence of an overdensity of gas in the 
vicinity of the quasar. 



the quasar. The only difference between the gas located either close 
to the quasar or far away from it is the increased photo-ionization 
rate in the vicinity of the quasar. In that case, far away from the 
quasar, the median optical depth corrected for redshift evolution, 
Tmedian, should be a constant we can call tq. In the vicinity of a 
QSO, due to the emission of ionizing photons, Tmedian is no more 
a constant (see Fig. 7) and decreases when r decreases as: 



Tmedian — Tq 



(9) 



where w was defined and derived in the previous Section (see 
Fig. 8). 

The observed median optical depth Tmedian is compared to 
To/(l + w) in Fig 9. It is apparent that the two curves do not agree 
and that the observed median optical depth is much larger than what 
would be expected in the case the density field of the IGM remains 
imperturbed when approaching the quasar. This means that on an 
average, there is an overdensity of gas compared to the IGM in the 
vicinity of the quasar. 

If an overdensity, is present around the quasar, then the 
median optical depth should be: 



Tmedian — To 



(pW/po) 



[2-0.7(7-1)1 



(I+W) 



(10) 



The 7 exponent has been shown by Schaye et al. (2000) to be in 
the range 7 = [1-1.5], thus it is reasonable to assume 7 = 1 for this 
work. 

Using Eq. 10 we can derive the actual overdensity. It is given 
versus the distance to the quasar in Fig. 10 (filled circles). The gas 
density close to the quasar is significantly higher than the mean 
density in the IGM at least within the first 10/i~^Mpc or so. 
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Figure 10. The mean density profile around quasars with emission redshift 
in the range 4 < Zcm < 'l-S as derived from this work is shown as filled 
circles. Open diamonds indicate the lower limit on the density profile de- 
rived assuming that the redshift from colimm 3 Table 1 is systematically 
vmderestimated by 1500 km s~^ and that the ionizing rate in the IGM is 
that measured by Becker et al. (2006). 



6 DISCUSSION AND CONCLUSIONS 

In this paper we have used the method presented by Rollinde et al. 
(2005) to probe the density structure around quasars. In the vicinity 
of the quasar and in comparison with the situation far away from it, 
the ionization factor of the gas is increased by the emission of ion- 
izing photons by the quasar and is decreased by the presence of an 
overdensity. If the redshift evolution of the IGM photo-ionization 
rate and that of the median optical depth in the IGM are both con- 
strained by the observed evolution of the median optical depth in 
the Ly-Q forest far from the quasar, then it is possible to derive the 
density structure in the vicinity of the quasar. We have followed 
McDonald & Miralda-Escude (2001) to derive the mean ionization 
rate in the IGM at 2; > 4. Our results are lower, by a factor of about 
1.5, than the new results by Becker et al. (2006). We then have used 
the redshift evolution of the pixel optical depth PDF to estimate the 
redshift evolution of the median optical depth of the IGM. Using 
these results, we have found that quasars are surrounded by signifi- 
cant overdensities on scales up to about 10 ft~^Mpc. The overden- 
sity is of the order of a factor of two at 10 ft~^Mpc but is larger 
than five within the first megaparsec. If true, we can estimate that 
the mass surrounding the quasars within 1 Mpc at this redshift is of 
the order of lO'^'* M© corresponding to the mass of a big cluster of 
galaxies. 

An overdensity can be artificially derived from the above analy- 
sis if (i) the photo-ionization rate is underestimated and/or (ii) the 
redshift of the quasars are systematically underestimated for any 
reason. To check the robustness of our result we can try to obtain a 
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lower limit on the overdensity surrounding the quasars. For this we 
have arbitrarily increased the photo-ionization rate that we used by 
a factor of 1.5 so that it is equal to the value obtained by Becker et 
al. (2006) at the same redshift. In addition, we have increased all 
quasar emission rcdshifts by 1500 km s^^ to account for a possi- 
ble systematic underestimate of the emission redshift when using 
the Ly-a emission line (note however that in our treatment we al- 
ready increase the emission redshift from the Ly-a emission line 
by an amount taken randomly between and 1500 km s~^). The 
result is shown in Fig. 10 (open diamonds). It can be seen that the 
overdensity is still present although smaller as expected. 
If our interpretation of the proximity effect is correct, we should 
expect a correlation between the strength of the observed prox- 
imity effect and the intrinsic luminosity of the QSOs. To test 
this effect, we have considered two subsets of our quasar sam- 
ple, one containing the 15 QSOs with highest intrinsic luminosi- 
ties (6.14 X 10+^^ < L < 1.22 X 10+^^ h'^erg s"^Hz-i) and 
the other containing the 17 QSOs with lowest intrinsic luminosities 
(1.85 X 10+^^ <L< 3.74 x 10+^^ ft-^erg s-^Hz"!). We have 
avoided quasars with intermediate luminosities because they may 
dilute the possible result. In Fig. 11, we plot for both subsets the 
median optical depth versus the distance to the quasar. It is apparent 
that the proximity effect is less pronounced for the low-luminosity 
sample compared to the high-luminosity sample as expected. The 
Proximity effect is correlated with QSO luminosities. For illustra- 
tion, we plot also T()/(1 + cj) for both subsets on the figure. The 
corresponding overdensities around the quasars in the two subsets 
are shown in Fig. 12. The overdensity is correlated with luminosity. 
Brighter quasars show higher overdensities. Some caution should 
be applied when interpreting this result however as we have used a 
quite simple model to correct the observed proximity effect. 

At redshift z ^ 2, Rollinde et al. (2005) claimed tenta- 
tive detection of overdensities of about a factor of two on scales 
~5/i~^ Mpc. The result was only marginal because the statistics 
was small. With a sample of 45 quasars at z > 4, we have demon- 
strated that overdensities exist. This result strongly supports the 
idea that quasars at high redshift are located in regions of high 
overdensities probably flagging the places where massive clusters 
of galaxies will form. 

It has been known for a long time that quasars are associated 
with enhancements in the distribution of galaxies (Bahcall, Schmidt 
& Gunn 1961) and that at low redshift (z < 0.4) they are asso- 
ciated with moderate groups of galaxies (e.g.. Fisher et al. 1996). 
However, little is known at high redshift. Our method yields a mean 
one-dimensional profile of the density distribution on large scales 
when other methods study the correlation of the quasar with sur- 
rounding objects (see Kauffmann & Haehnelt 2002). This result 
adds to the growing evidence that high-z QSOs seem to reside in 
dense and probably highly biased regions (Djorgovski 1999, Djor- 
govski et al. 1999). Observational efforts should be done to obtain 
deep images of the fields around bright high-z quasars to search for 
the presence of any concentration of objects aroirnd the quasar. 

Acknowledgements: RG is supported by a grant from the brazil- 
ian gouvernment (CAPES/MFC). AA is supported by a PhD grant 
from the Ministry of Science, Research &Technology of Iran. PPJ 
and RS gratefully acknowledge support from the Indo-French Cen- 
tre for the Promotion of Advanced Research (Centre Franco-Indien 
pour la Promotion de la Recherche Avancee) under contract No. 
3004-3. SGD is supported by the NSF grant AST-0407448, and 
the Ajax Foimdation. Cataloguing of DPOSS and discovery of 




0.01 



High luminosity opticai depth i — B- 
Low iuminosity opticai deptin — 
TQy(1+(D) liigli iuminosity r--^- 
Tq / (1+0)) iow iuminosity 



10 20 30 40 50 60 70 

Distance (Mpc.h ^ 

Figure 11. The median optical depth in versus the distance from the cloud 
to the quasar for two subsets of our QSO sample: one including the 15 QSOs 
of highest luminosities (filled squares) and the other including the 17 QSOs 
of lowest luminosities (open triangles). It is apparent that the proximity 
effect is correlated with luminosity as expected. The curves giving the factor 
(1 +u}) are drawn for each subset for illustration (filled circles and open 
diamonds respectively. 
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